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ABSTRACT 

A  method  is  presented  for  the  solution  of  second 
order  electromagnetic  fields  in  a  bounded  plasma  which  is 
excited  by  two  high  frequency  electromagnetic  waves.  In 
particular,  the  problem  considered  is  the  optimization  of 
second  order  fields  in  a  uniform  plasma  confined  to  a 
circular  guide.  Consideration  is  also  given  to  the  amount 
of  ion  heating  that  may  be  expected  through  the  wave  mixing 
process  when  a  static  magnetic  field  is  applied  to  the 
plasma.  Two  underlying  assumptions  in  the  analysis  are: 

(1)  the  frequencies  of  the  incident  waves  are  sufficiently 
high  that  they  achieve  a  thorough,  uniform  penetration  of 
the  plasma;  (2)  the  thermal  velocities  of  the  species 
comprising  the  plasma  are  much  less  than  the  applicable 
phase  velocities  of  any  disturbances;  the  latter  assumption 
justifying  a  cold  plasma  model. 
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CHAPTER  1 


I  NTRODUCT I  ON 


1.1  General 

In  recent  years,  the  problems  associated  with  the 
nonlinear  interaction  of  high  frequency  electromagnetic 
waves  in  a  plasma  have  drawn  much  attention.  It  has  been 
found  in  a  number  of  theoretical  papers  that  it  is  possible 
to  achieve  a  resonant  mixing  of  two  waves  of  frequencies  to  ^ 
and  002  propagating  in  an  unmagnetized  plasma  when  the  dif¬ 
ference  frequency,  to  ^  -  to^  ,  is  approximately  equal  to  the 
e  1  e  c  t  r  o  n  -  p  1  a  s  ma  frequency  tOpe  =  [^TTe2nQ/me]  .  The  mixing 
process  can  result  in  appreciable  second  order  fields. 

In  1964,  Salat  and  Schluter  (1),  described  a 
diagnostic  procedure  whereby  plasma  density  could  be 
determined  by  measuring  the  intensity  of  light  scattered 
by  a  plasma  excited  by  a  single  laser  beam.  The  scattered 
light  rays  from  different  points  in  the  scattering  volume 
were  generally  out  of  phase  except  in  a  particular  scat¬ 
tering  direction.  By  measuring  the  angle  between  the 
inc'ident  signal  and  the  direction  in  which  light  was 
scattered  in  phase  at  the  harmonic  frequency  2c0|  ,  the 

plasma  density  could  be  determined.  A  major  disadvantage, 
however,  was  that  the  amount  of  energy  scattered  was  ex¬ 
tremely  small  in  relation  to  the  source  energy. 

A  modification  of  the  above  procedure,  using  two 
laser  beams,  results  in  a  much  more  favourable  scattered 


2 


to  incident  energy  ratio.  Such  a  method  was  considered  by 
Kroll  et  al  (2).  By  mixing  two  laser  beams,  it  was  found 
that  considerable  enhancement  of  the  scattering  cross-section 
was  achieved  when  the  difference  frequency  co  ]  -  002  -  Wpe  , 

giving  a  significant  amplitude  to  the  scattered  wave.  A 
comparison  was  drawn  for  a  specific  example  showing  that 
the  level  of  the  scattered  signal  using  a  single  laser,  as 
compared  to  using  two  lasers,  was  -  10  16.  A  more  detailed 

account  of  this  procedure  is  given  by  Shkarofsky  (3). 

A  further  modification  has  been  suggested  by 
introducing  a  third  laser  beam.  A  disadvantage  of  the 
above  methods  is  that  the  detectors  must  look  almost  direct¬ 
ly  into  the  incident  beams  because  the  angle  between  the 
scattered  light  and  the  laser  beams  is  quite  small.  Salat 
(4)  has  shown  that  scattering  a  third  beam  off  the  fluctu¬ 
ations  due  to  the  first  two  can  lessen  this  problem 
considerably,  however  the  orientation  of  the  third  beam  is 
subject  to  high  tolerances. 

The  above  diagnostic  procedures  have  involved  the 
use  of  laser  beams  and  treated  the  plasma  as  an  infinite 
medium.  It  has  recently  been  found  by  Stern  and  Tzoar  (5), 
that  two  microwave  signals  will  stimulate  oscillations  in  a 
plasma  at  electron  densities  at  which  the  plasma  can  support 
resonant  modes  (for  a  uniform  plasma,  this  is  the  electron- 
plasma  frequency).  With  this  type  of  problem,  where  the 
incident  waves  mix  throughout  the  plasma  and  are  a  function 
of  the  confining  structure,  the  plasma  must  be  considered 
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as  finite  to  achieve  a  correct  interpretation  of  experimen¬ 
tal  results.  In  this  thesis,  a  theoretical  approach  to  the 
problem  of  wave  interaction  in  a  finite  plasma  is  presented. 
Initially,  a  treatment  of  an  unmagnetized  plasma  will  be 
given,  followed  by  a  similar  problem  when  a  static  magnetic 
field  is  applied  in  the  direction  of  propagation  of  the 
incident  waves.  With  respect  to  the  latter  case,  it  has 
been  suggested  by  James  and  Thompson  (6),  that  a  possibility 
exists  for  significant  ion  heating  through  a  wave  mixing 
process.  However,  their  theory  was  based  on  an  infinite 
plasma  excited  by  plane  waves.  To  obtain  a  more  accurate 
picture  of  what  is  achievable,  the  plasma  must  be  considered 
as  finite,  since  the  degree  of  coupling  will  have  a  depen¬ 
dence  on  the  shape  of  the  confining  structure.  Hence  a 
major  consideration  for  the  magnetized  plasma  will  be  an 
investigation  into  the  degree  of  ion  heating  that  might  be 
expected  in  a  finite  structure. 

1.2  Plasma  Model 

In  this  thesis  the  cold  plasma  model  will  be 
used.  It  will  be  assumed  that  each  charge  species  (ions 
and  electrons)  obeys  an  equation  of  continuity: 

3_n_a  +  V*(nav.a)  =  0  ...(1.1) 

3  t 


and  an  equation  of  motion: 


>- 
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9  ,  . 

—=za  +  (  v  •  V  )  v 

TTF  —a  —a 


ia  (  E  +  ZaJL 

ma  V  —  c 


+  ^(collision  terms) ^  ...(1.2) 


a 


Terms  involving  collisions  will  be  considered  later,  depend¬ 
ing  on  the  particular  problem  being  investigated. 

In  an  unmagnetized  plasma,  thermal  effects  can  be 
neglected  provided  (7)  the  thermal  velocities,  v^  ,  of  the 
ions  and  electrons  are  much  less  than  the  phase  velocities 
of  any  disturbances. 


v  0  a 


...(1.3) 


where  oo  and  k 
di sturbance. 
plasma  is 


are  the  frequency  and  wavenumber  of  the 
The  equivalent  condition  for  a  magnetized 


...0.4) 


where  ft ^  is  the  cyclotron  frequency  of  the  a-species. 
Further,  it  has  been  verified  (8),  when  (1.3)  and  (1.4) 
are  satisfied,  that  to  second  order  the  density  and  current 
fluctuations  produced  in  a  plasma  by  external  waves  are  the 
same,  whether  derived  from  kinetic  theory  or  cold  plasma 
theory.  It  follows  that  the  resulting  fields  must  neces¬ 
sarily  be  the  same  from  either  approach.  Since  cold  plasma 
theory  is  mathematically  less  involved  than  kinetic  theory, 


■ 
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this  approach  has  been  used. 


1.3  General  Equations 


Assume  two  electromagnetic  waves  of  frequencies 


cu  i  and  o)2  are  incident  on  a  uniform  plasma  confined  to  a 
cylindrical  waveguide  of  circular  cross-section.  In 
addition,  assume  that  to  ]  and  002  are  sufficiently  high  that 
shielding  of  these  waves  is  negligible  to  first  order  - 
there  is  essentially  an  impedance  match  between  free  space 
and  the  plasma.  Hence,  no  reflection  of  the  incident  waves 
occurs  and  the  waves  are  not  attenuated  in  the  plasma. 

Under  this  assumption,  the  first  order  fields  are  the 
solutions  of  Maxwell's  equations  in  an  empty  waveguide. 


To  consider  a  particular  problem,  assume  that  the 


incident  waves  are  transverse  electric,  TEQm  and  TEQn, 
modes.  Hence,  from  Maxwell's  equations: 


V  •  E 


4  TT  p 


...(1.5) 


V  •  B 


0 


V  x  E 


I  3  B 
c  3  t 


...(1.7) 


V  x  B 


iLl  i  +  1  H 

c  c  3 1 


and  using  cylindrical  co-ordinates  (  r  ,  tf>  ,  2  )  ,  the  solutions 


for  the  incident  fields  are: 
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(a)  TE 


0  m 


B  o  (  k  i  |r)cos(w^t  +  k2|Z  +  ap 

E^J 1(k] ]r)sin(w]t+k2|z+a|) 

B  *J  !  (  k  1  |r)sin(w1t  +  k2|z  +  a]) 


(1.9) 


where 


1  1 


m 


...(1.10) 


h 

with  s  |  m  being  the  m L  root  of  Jj(s)  =  0,  and  rQ  the  radius 
of  the  waveguide.  Also 


Bz 


1  1 


+  k 


to 


2 

2  1 


k"c  e| ,  bJ. 


k21° 


CO 


..(1.11) 

..(1.12) 


(b)  TE 


0  n 


B  z 

El 

B2 

r 


BzJo(k12r^cos^w2t+k22Z+a2' 

E  ?  J  i  (k12r)sin(w2t  +  k22z  +  a2) 
ErJl  (k12r)sin(oj2t  +  k22z  +  a2) 


(1.13) 


where 


1  2 


...(1.14) 


w 


ith  s1n  being  the  nth  root  of  J|(s)  =  0.  Also 


2  2 
k  1  2  +  k2  2 


..(1.15) 


~B2z 


1  2C  E  ? 


bI 


CO. 


k22C  E| 


CO. 


(1.16) 


The  fields  are  now  expressed  in  exponential  form: 
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B  1 

V 

B  1  “ 

z 

z 

z 

E  \ 

B  1 

i  (ajj  t  +  k2  j  z) 

E1 

% 

B  1 

-i (wi t+k,, z) 

+  e  1  z  1 

✓N  -  J, 

E  4> 

B  1  " 

r 

r 

r 

where 


...0.17) 


H  -  I  BzJ0(k)  1  r)e  1011 

B  \  =  I  B‘j)(k))r)e'C<1 


B2 

A  r\ 

B2 

a  .V 

B2 

z 

z 

z 

E* 

B2 

1 

<p 

B2 

+  e "  '  (w2t  +  k22z) 

E* 

B2* 

r 

r 

r 

where 

*1  ■  i  B2j0(k,2r)eia2 

2 

E2  =  1  E2J  (k  r)e'a2 

b  4)  '  1  ^k12r;e 
2  i 

B2  =  1  B2 J ]  ( k ! 2  r ) e ' a2 


(1.18) 


...(1.19) 


\ 

(1.20) 


CHAPTER  2 


NONLINEAR  INTERACTION  OF  WAVES  IN 


AN  UNMAGNETIZED  PLASMA 

2.1  Introduction 

It  is  assumed  in  this  chapter  that  a  TEQm  and  a 
TEq0  wave  of  the  type  discussed  in  Art.  1.3  propagate  in 
an  unmagnetized  plasma  confined  to  a  cylindrical  waveguide. 
Since  E_  and  B_  are  unknown,  equations  (1.1)  and  (1.2)  must 
be  solved  in  a  self-consistent  manner  with  Maxwell's 
Equations,  hence  an  extremely  complicated  set  of  nonlinear, 
coupled  partial  differential  equations  must  be  dealt  with. 
Due  to  the  n o n 1  i n e a r  i  t i e s  ,  superposition  principles  do  not 
hold  (the  plasma  mixes  the  incoming  waves  through  the  terms 
(v_*V)_v  and  j_xB_)  ,  hence  sum  and  difference  frequency  waves 
are  generated.  Under  certain  conditions,  significant 
oscillations  are  stimulated  at  the  difference  frequency, 
w  ]  -  0J2  •  To  investigate  the  nature  of  these  conditions  and 

the  size  of  the  oscillations,  the  equations  describing  the 
plasma  are  linearized  by  assuming  an  expansion  to  second 
order  of  the  form: 


n 

~no~ 

^l' 

V1  ‘ 

V 

0 

*1 

E1 

V? 

E1* 

E. 

- 

0 

+  e' (“lt+k21z) 

+  e_i (“lt+k21z) 

B_ 

0 

B  1 

B  1  * 

- 

.. 

—  _ 
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n  2 

r  -a-  1 

-  /v 

n  2 

i  ( co 2  t  +  k2  2z  ) 

—2 

E_2 

+  e-i  (a,2t  +  k22z) 

Ju 

—2 

E_2* 

B2 

B  2  * 

n  3 

*  JU  • 

n  3 

i  ( cot  +  kz  ) 

+  e 

—3 

L3 

-  i  ( co  t  +  kz  ) 

+  e 

4 

^  3  * 

I3 

B  3  * 

where 


co  =  co  ^  -  co  2 


k  =  k 


2  1 


-  k 


..(2.1) 
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.  .  .  (2.2) 

.  .  .  (2.3) 


In  equation  (2.1),  the  functional  dependence  of 
n,  v,  E_,  and  B_  on  r  is  understood.  First  order  terms  are 
denoted  by  1  or  2  -  i.e:  terms  which  vary  in  time  at 
frequency  co  |  o  r  •  Second  order  terms  are  denoted  by  3* 


2.2  Collision  Te  rms 

The  primary  interest  in  this  chapter  is  to  stimu¬ 
late  oscillations  near  the  e 1 e c t r on - p 1  a s ma  frequency,  hence 
it  is  assumed  that  the  difference  frequency,  co ,  is 
sufficiently  high  that  the  motion  of  the  ions  can  be 
neglected.  Equation  (1.2)  then  reduces  to  an  equation  of 
motion  for  the  electrons  of  the  type  considered  by  Salat  (3) 


3- 

•  • 


(2.4) 


-10- 


3v_e  +  (^ev))ie  =  0_e(I  +  YJi)  “  vei^_e 
3 1  me  c 


where  ve j  is  the  electron-ion  collision  frequency. 


2.3  First  Order  Fields  and  Velocities 


Since  it  has  been  assumed  that  the  incident  fields 
are  unaffected  by  the  presence  of  the  plasma,  comparing 
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Equating  terms  of  dependence  e  z  gives: 
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2.4  Derivation  of  Second  Order  Currents 


Second  order  terms  are  considered  by  equating 
terms  of  dependence  e'(wt+kz)  when  equation  (2.1)  is 
substituted  into  equation  (2.4).  The  result  is; 

iU-ivei)v_3  +  (  *  v  t )  v.2  "  '  ^  if_2  *  —1  ^  —2 

+  ( ^_2  '  V  t  ^  —  1  +  '  1  *-P— 1 


where 


me  c 


vf  B_2* 


— z  k  2  1 


— zk2  2 


...(2.11) 


V  denotes  the  component  of  the  V-operator  which  is 
perpendicular  to  u_^. 

Making  use  of  equations  (2.5)  to  (2.10),  equation 
(2.11)  is  now  considered  term  by  term. 
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\  E  1  •  V  t )  E_2  * 

% 

Lme, 

(O  |  to  2 

me 

where  R  =  -  1  E  |,E  | 

to  I  a)  2  r 


b  ]  (  k_2  *  v_i  )  v_2 


1  qe 


2 


k22 

(0  |  (O2 


0 


c]  (v|*Vt)v1  = 

^  e 

2  1  U2*‘Vt)E^  = 

\?i 

1 _ 

_me_ 

a)  |  0)2 

1 

3 

CD 

1 

...(2.11a) 

...(2.11b) 

...(2.11c) 


c]  (v|*Vt)v1 


•  •  • 


(2.  1  Id) 


Jv< ,V* J v) 


-12- 


d  J  ^ili  'Li  )-l 


1  2 


m 


k2  1  (u.z  'I2*)!1  =  0  ...  (2  .  1  1  e) 


e ]  q  e  (^jxB_2“)  =  " 


mec 


m. 


1  (u  eTb2* 
_  —  r  <p  z 


0)  |  c 


-  u  E  !  B  2  *  ) 
— z  q  r 


q  e  (V2xi1}  =  1 


-  n  2 


mec 


m. 


1  <^r§z 


u)2  c 


A  kU  A  i  v 

0zE|  Br> 


Hence 


q  e  (v_lxi2"  +  v_2x  i1  )  = 

<i)| 

1 

2 

A  A 

S  u  r  +  Qu  z 

mec 

m  e 

.  .  (2 .  1  1 f ) 


where 


S  =  -  i 


p  1  S  2  *  _ 

E  4>Bz  e  <j)  8  z 


0) 


0). 


Q  = 


e i b ^ -  e|"b‘ 


.  .  .  (2 .  1  lg) 


. . . (2 . 1 1 h) 


0) 


0)  ■ 


Substituting  equations  (2.11a)  to  ( 2 .  1  1 h )  into  equation  (2.11) 
gives: 

i  (oj  -  i  ve  j  )  v 


v3  =  qeE_3  +  jq  e  ]  2  ^iLr 


me  |me  ! 


+  It 

me 


2? 


where 


P  =  S  -  2$ 


Q  u  z 

..(2.12) 

..(2.13) 


It  can  be  readily  shown  that: 

P  =  "  1  _ |_EiE6e '  (ai  a2^  3 [ J i  ( k i  i r ) J i  (k1?r)] 


co  |  0)2 


8  r 


(2.14) 


and 


A  --  I  —  O  I  (  Ot  I  “  Ot  Q  )  /  \  /  \ 

Q  =  -  j.  _ L_EJE<j>e  kJ  )  (k)2r)  J  !  (k,  ,  r) 

4  co  ]  o)2 


(2.15) 


13 


Note, 


and  ax 


or,  in 


where 


and 


and 


from  equations  (2  . ] b)  and  (2.15): 


/\ 


=  i  kP 
3  r 


...(2.16) 


Reducing  equation  (2.12)  to  its  radial,  angular 
ial  components  gives: 


i  ( to 
i  (  <0 
i  ( a) 


v 


ei)v3r  =  ^eEr 


m, 


ivei)v3* 


-  i  v  .  )  v 


i eEl 


e  i  3  z  M  e 


m, 


q  E3  + 
M  e  z 


m. 


m. 


2P 


m. 


2$ 


(2.17) 


J 


terms  of  currents: 


-1  3  r 

j34> 
J*  3z 


a  E  3  +  b  P 
^Tri  b  tt  i 


a 

b  it  i 

a  E 3  +  b  (J 


(2.18) 


b  TT  i  An  i 


j  0  =  n  q  v  -> 
J-3  oM  e— 3 


(2.19) 


a  —  to 


pe 


(0-  i  ve  j 


b  =  - 


2 

wpe 


m. 


( w- i vP  :  ) 


(2.20) 


P  =  |  e  |  P 

$  =  lei  6 


(2.21  ) 


-  1 4  - 

Note  that  equation  (2.19)  is  valid  because  n  j  and 
r>2  are  both  zero;  the  plasma  density  is  not  perturbed  to 
first  order.  This  follows  from  the  equation  of  continuity 
(1.1)  for  electrons. 


a n e  +  V-  (neve) 

=  0 

...(2.22) 

at 

Substituting  equation  (2. 

1  )  into 

equation  (2.22) 

and  equa  t  i 

terms  of  dependence  e '  ^ w 1 

t  +  k2  j  z ) 

gives: 

i  (jo  i  n  i  +  nQ  (  i  k_2  i  *v_)  )  +  nQ(Vt»vj)  =  0  ...(2.23) 

Since  v_|  is  in  the  ^-direction  only,  and  is  not  a  function 
of  4>  or  z,  then  the  last  two  terms  of  equation  (2.23)  are 
zero.  Hence  n^  =  0.  The  same  is  true  for  n2* 


2.5  Differential  Equations  for  Second  Order  Fields 

From  Maxwell's  Equations,  the  following  relation 
ships  can  be  derived  (see  Appendix  1). 


3_ 
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F  3 
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=  i  4  it  (jo 


-134) 


I  i_  rlLz 

r  a  r  V  a  r 
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a»\2-k2 


e  i 


2  . 


2  . 


=  ?  4  IT  j  -  k  j  3  z  +  fw\"j3z  +  ik  1  9  (  r  j  3  r  ) 

oo  [_  \C  )  r  a  r 


(2.24) 
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Substituting  equations  (2.18) 


v  w  3rlr  3r  / 


b_ 

(0 


i_/i  i_(r"p)  + 

3  r  r  3  r  ] 


1_/1  3_(rE|)\  + 

3  r  \  r  3  r 


1  -_a_ 

0) 


I  <L(  rllz 

r  3r\  3  r  ; 

=  b_ 

0) 


00  \2-k2  w  1  - a_ 

00 


oo  \  2  -  k  2  ]  Q  +  ikl 


into  equations 


oo  / 


oo_\  2  P  +  i  k  3  Q 
c  j  3  r 


(  2  .  2  A  )  gives: 

i  ka  3_eJ 

to  3  r 

... (2.25) 


...(2.26) 


Ez  _  iki  I  3_(  rE3) 

oo  r  3  r 


9_(rP) 
3  r 


.  .  .  (2.27) 


Consider  equation  (2.26).  There  is  neither  a 
driving  function  to  excite,  nor  coupling  of  the  other  field 
components,  E3  and  E3,  to,  E|.  The  only  possible  solution 
is  the  homogenous  solution,  however,  since  initially  there 
are  no  second  order  fields,  then  =  0.  Hence,  using 
equation  (2.16),  the  set  of  equations  to  be  solved  reduces 
t  o  : 

i  ka_  3  E  3 
to  3  r 

...(2.28) 

]_  3_/r3_E|\  +  //w)  2-k2W  1  -a\  E3  -  i  ka_  J_  3_(  r  E  j! ) 
r  3  r  \  3  r  I  \\c  I  l\  ul  to  r  3  r 

... (2.29) 


=  b_ 

3 

/r  3 Q  \  + 

r 

1  o' 

CM 
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CM 

3 1 

to 

r 

3  r 

\  3r  J 

\C  / 

(  1  9_(  rE  r  )]  + 

\  <o/3r\r  3r  /  \c 


=  b_ 

to 


-  k 


2\  s 


E  3  - 
c  r 


' 


1  6 


2.6  Solution  of  Second  Order  Field  Equations 


The  approach  used  to  solve  for  the  fields,  and 

E|,  involves  expanding  E|  into  a  complete  set  of  orthonormal 
functions  which  are  consistent  with  the  boundary  condition, 
E|(r=r0)  =  0.  Taking  into  account  the  form  of  equation 
(2.29),  a  Fo u r  i  e r - B e s s e  1  expansion  of  zeroeth  order  Bessel 
functions  appears  appropriate.  Hence  it  is  assumed  that: 


E  |  (  r )  =  l  A  j  J  0  (  r  £  ,  ) 
i  =  1 

where 

Jo(ro^i )  =  0 

o  r 


...(2.30) 
. .  .  (2.31  ) 


*i 


=  s 


o  i 


r 


.  .  .  (2.32) 


o 

where  s0j  is  the  i ^ ^  zero  of  a  zeroeth  order  Bessel  function. 
The  problem  remaining  is  to  determine  the  Fourier-Bessel 
coefficientsAj. 

Substituting  equation  (2.30)  into  equation  (2.29) 


gives  : 

oo 

i  ka_  J_  9_(  r  E  j* )  =  £  A;r;J0(r£j)  -  R2(r) 

oo  r  8  r  i  =  1 

where 


=  b_ 

1  3  ( 

r  9  Q. 

+ 

i 

CM 

1 

CSI 

3 

Q 

w 

r  8  r 

\  dr  j 

\yj 

...(2.33) 


...  (2.3*0 

.  .  .  (2.35) 


Integrating  equation  (2.33)  from  0  to  r  gives: 


' 


. 


1  7 


00  a  •  r  • 

ii 


00 


1  =  1  M 


iklEr  =  I  -^-L  j  l  ( re  I  >  -  i  Jr  '  R2  (  r  '  )dr  '  ...  (2.36) 


0 


Note,  from  equation  (2.30)  : 


3Ej 

3  r 


l  A  i  s  I  J  1  (  re  i  ) 

i  =  1 


.  .  .  (2.37) 


Substituting  equations  (2.30)  and  (2.37)  into  equation 
(2.28)  gives: 


A  i  A  j 


1  £ 

i  =  1  5 


J  1  ( r  £  j )  =  i ka  R !  ( r)  +  (\ -a\ 9R2 ( r> 


,ai\| 

00 
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-  k‘ 


oo  3  r 
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where 


and 
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P 
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- 
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1  -a_ 

oo 


_  2  .  2 
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oo)2  -  Si 

vC 


-  k' 


...(2.40) 

Multiplying  equation  (2.38)  by  r  and  differentiating  with 
respect  to  r  gives: 

...(2.41) 


£  A  j  A  j  J  o  (  r  £  j  )  =  S  (  r  ) 

i  =  1 


where 


S  (  r  )  = 


1  -a_ 

oo 


-  k‘ 


R  2  (  r ) 

...(2.42) 
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Using  the  ort hogana 1 i ty  property: 
ro 

/rJ0(r£j)J0(r£j)dr  =  0 


0 


-  Lo  Jl<ro5j) 
2 


i  *  J 


'  =  J 


It  may  readily  be  shown  that: 


A 


2  1 


roAiJ?<roSi)  0 


JrJ0(r5|)S(r)dr 


.  .  .  (2.43) 


Hence,  from  equation  (2.30): 

2  00  J n  ( r?  i  )  r° 

Ez(r)  =  “2  I  ~  72~7 - TT  /rJ0(r?i )S(r)dr 

roi  =  l  Ai J1  (ro^i  }  0 


.  .  .  (2 . 44) 


and,  from  equation  (2.36) 
E  3  (  r)  =  w  _J_ 


a  i  k 


2  “  r ; J i (rg;)  o 


“  I 

ro  i=1  ^iAiJ?(ro^i)  0 


JrJ0(r^j)S(r)dr 


I  |r 1 R2  (  r  1  ) dr 
r  0 


.  .  .  (2.45) 


2.7  Optimization  of  Second  Order  Fields 


As  mentioned  earlier,  a  resonance  effect  near  the 
e 1 ec t ron - p 1  a s ma  frequency  has  been  predicted  theoretically 
for  an  infinite  plasma.  It  will  now  be  shown  how  the  reson¬ 
ance  phenonema  can  be  accounted  for  in  a  bounded  plasma. 

Consider  the  expression  for  A;  -  equation  (2.40). 
Both  E|  and  are  proportional  to  1/(1— a/oi)  .  From  equation 
(2.20): 


£ 

a) 


“Pe  (\  +  i 

co2  +  v|  i  \  co  J 


...(2.46) 
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Assuming  u>  >>  vej 


£ 

0) 


. . . (2.47) 


I  f  to  =  Upe  ,  then  1  -  a_  =  -  i  v  e  i  .  Since  E|  and  are  then 

w  wpe 

proportional  to  Wpe/vej,  it  can  be  seen  that  the  possibility 
exists  for  a  significant  peak.  In  general,  a  maximum  will 
occur  when  : 


to  I  =  ^pe 


e  i 


co 


pe 


. . . (2.48) 


Fig.  2.1 


co  /  v  .  , 
pe  ei  ’ 


shows  a  graph  of  |  1  -  a/co| 

in  the  region  to  'v  0(10^). 


for  various  values  of 


2.8  Sensitivity  of  Field  Solutions  Near  the  Optimum 
Point 


Since  the  factor  1/(1  -  a/co)  is  the  dominant  term 

affecting  the  magnitude  of  the  second  order  fields  near  the 
resonance  point,  a  measure  of  its  variation  with  the  dif¬ 
ference  frequency  to,  and  the  plasma  density  nQ,  will  give 
a  good  indication  of  the  overall  sensitivity  of  the  fields 
to  fluctuations  in  these  parameters. 

When  w p e  <  10vej  the  resonance  peak  is  not  too 
sharp  and  sensitivity  will  not  be  a  serious  problem,  how¬ 
ever,  when  tOp6  >  1  0  v  e  j  ,  the  peak  can  become  very  sharp.  In 

general,  the  optimum  condition  can  be  stated  as:  | X+ i Y  |  be 
a  minimum,  where: 


20 


-21- 


X  =  Re(l-a/to),  V  =  Imfl-a/w) 

This  condition  is  satisfied  when  equation  (2.48)  holds.  If 
Upg  >  10vej,  then  the  optimizing  condition  can  be  stated  as 

X  ( to ,  n  0 )  =  0  ...(2.49) 


If  the  optimum  condition  is  not  to  be  departed  from  too 
greatly,  then  the  following  must  hold: 


|  X  ((jo+6co  ,  n0+6n0 


3X 


3  to 


6  to  + 


3  X 


3  rfc> 


6  n , 


...(2.50) 


If  the  inequality  is  considered  along  each  axis,  then 


I  | 


3_X 

3  uo 


...(2.51) 


6  n 


3X. 

3 1\> 


...(2.52) 


Since  3_X  -  2  /^p  e  \  J_,  and  3_X  -  -  ]_  ( ^pe  \  ,  and  since  to-tope 
3 oo  \  to  )  to  3nQ  no\  w  / 

equations  (2.51)  and  (2.52)  become: 


6 oj  I  ^  ve 


.  .  .  (2.53) 


6  n, 


jc  1  ve  i  n 


2  to 


pe 


•  •  • 


(2.54) 
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2.9  Example  Problem  1  -  Complete  Solution  of  Second  Order 

Fields  in  an  Unmagnetized  Plasma 

A  complete  solution  of  the  average  value  of  second 
order  fields  as  a  function  of  the  difference  frequency,  co , 
is  now  given.  The  problem  was  solved  with  the  aid  of  an 
IBM  7040  computer  which  carried  the  field  expansions  defined 
by  equations  (2.AA)  and  (2.A5)  to  ten  terms.  The  tenth 
terms  were  'v  0(10  ^ )  times  the  first  terms  of  the  series 
and  the  sign  of  successive  terms  was  alternating  indicating 
strong  convergence  of  the  solutions. 

To  give  a  clear  indication  of  the  spectrum  of  the 
fields,  the  plasma  was  defined  as  follows: 

co  =  1.5X10^  rad/sec 


^  =  ]00 


The  wa  veg  u 

i  de  r  a  d  i 

us  was 

taken  as  2 

cm .  and  the 

firs 

t  i  n  c  i 

dent  field 

was  ass  umed  to 

be  a  T E o  , 

wave  with  a 

f  i  xe 

d 

frequency, 

co ,  ,  of 

9*  1  0  1  0 

r  a  d  /  s  e  c  . 

The  frequency 

o  f 

the 

second  inc 

i  d  e  n  t  f  i 

eld,  also  a  T  E  q 1 

wave,  was  var 

i  e  d 

from 

1 3* 1 0 1 0  ra 

d  /  s  e  c  to 

o 

X 

OO 

vO 

^  rad/sec 

,  below  which 

the 

field 

would  not  propagate  in  the  waveguide. 

The  magnitude  of  the  average  value  of  the  second 
order  fields  can  be  expressed  in  the  form: 


<E^>  =  ArU)|E;E2| 

<E3>  =  Az(u) I E‘E|| 


■(2.55) 
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Fig.  2.2  illustrates  the  variation  of  1 o g ^ q [ A r ( w ) ]  and 
1  og  ^  q  [  Az  ( oo )  ]  as  a  function  of  oi  giving  the  expected  reson¬ 
ances  when  |  w  |  =  ojpe  . 

2.10  Concluding  Remarks 

The  analysis  presented  in  this  chapter  has  verified 
for  a  special  case,  the  observations  of  Stern  and  Tzoar  (5); 
namely,  that  significant  oscillations  can  be  stimulated  in 
a  plasma  by  two  high  frequency  signals  whose  difference 
frequency  is  near  a  resonant  mode  of  the  plasma. 

The  use  of  wave  mixing  as  a  diagnostic  has  been 
supported  theoretically,  A  particular  advantage  of  the 
modes  selected  for  excitation  is  that  the  plasma  density 
is  not  perturbed  to  first  order.  This  will  not  be  true  in 
general  . 

Equations  (2.53)  and  (2.5M  indicate  a  rather  high 
degree  of  tuning  and  low  variation  in  density  would  have  to 
be  achieved  to  realize  an  optimum  condition,  if  the  electron- 
ion  collision  frequency  is  low.  However,  this  is  not  a 
necessity  when  the  procedure  is  used  as  a  diagnostic.  This 
is  apparent  from  Fig.  2.1.  The  line  shape  is  almost  iden¬ 
tical  in  the  lower  region  of  the  graph  for  all  values  of 
00pe/vei  >  100,  therefore  points  off  resonance  can  be  projec¬ 

ted  to  give  a  good  indication  of  where  the  actual  resonance 
will  lie. 


Ip  I 


CHAPTER  3 


NONLINEAR  INTERACTION  OF  WAVES  IN 


MAGNETIZED  PLASMA 

3.1  Introduction 

The  problem  which  will  now  be  examined  is  that  of 
nonlinear  interaction  of  waves  in  a  plasma  confined  to  a 
cylindrical  waveguide  with  a  static  magnetic  field,  B  , 
applied  axially  along  the  guide.  Again,  excitation  will  be 
by  high  frequency  T  E  q  m  and  T  E  n  o  modes.  In  contrast  to 


0  n 


Chapter  Two,  the  difference  frequency  will  be  significantly 
lower  than  go  ;  sufficiently  low  that  ion  motion  cannot  be 
neglected.  The  intent  here  is  to  investigate  what  degree 
of  ion  heating  may  be  expected  through  the  wave  mixing 


process 


The  magnetic  field  Bq  introduces  additiona 


characteristic  frequencies  defined  as  cyclotron  frequencies 


ft 


a 


-  9  a  B  o  , 


mac 


(3.1) 


is  the  cyclotron  frequency  of  the  a-speci es  . 

The  methods  used  to  solve  the  problem  are  similar 
to  those  employed  in  Chapter  Two,  therefore  many  of  the 
details  are  not  included.  The  analysis  is  begun  by  assuming 
an  expansion  similar  to  equation  (2.1),  however,  both  ions 
and  electrons  must  be  accounted  for. 
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It  is  desirable  that  the  second  order  fields  be 
at  the  cut-off  point  so  that  none  of  the  energy  available 
to  heat  the  ions  propagates  from  the  plasma.  This  means 
that  the  wavenumber  k  =  0.  If  this  condition  is  imposed 
with  the  additional  requirement  that  the  incident  waves 
propagate  in  the  waveguide,  it  is  possible  to  establish 
boundaries  which  restrict  the  range  of  oj  j  and  once  the 
modes  of  excitation  have  been  selected. 

If  k  =  0,  then  k^|  =  ^22'  F rom  equations  (1.10), 

(1.11),  (l.U)  and  (1.15) 


2 
r  ^ 


...(3.3) 


m 
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If  the  incident  fields  are  to  propagate  then 


Hence 


1 


o)2 


s  1  n 


I  m 


[s 


.  .  .  (3.4) 


u2  <  s_j_n 
W1  <  s  1  m 


[s 


1  m 


<  s 


(3-5) 


These  conditions  are  illustrated  graphically  in  Fig.  3.1. 
The  region  where  all  conditions  are  satisfied  is  shaded. 


Fig.  3*1  Regions  that  permit  second  order  cut-off 
waves  and  propagating  incident  waves. 
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3.2  Collision  Terms 


Primary  concern  in  this  chapter  is  the  effect  of 
wave  mixing  on  the  motion  of  the  ions.  The  difference  fre¬ 
quency  w  is  sufficiently  low  that  ion  motion  is  significant, 
hence  an  equation  of  motion  for  both  ions  and  electrons  must 
be  considered. 

-  ve  j  (  v_e  -  V|)  ...(3.6) 

-  \^)  ...  (3  .7) 


9  v_i 
Tt 


+  (  V  ;  •  V  )  V 


Vi  x~ 


B  \  -h  v  P  ;  (  v 


e  i  '  __e 


^e 

3  t 


+  (ve-V)ve  =  q 


e 


ve 


where 


m  r  v  e 


(3.8) 


3.3  First  Order  Fields  and  Velocities 

If  the  additional  assumption,  co  |  ^  >>  ^  e  >  's  made, 

then  the  equations  for  the  first  order  fields  are  identical 
to  equations  (2.5)  to  (2.8). 

The  equations  for  first  order  velocities  are 
similar  to  equations  (2.9)  and  (2.10),  but  both  ions  and 
electrons  must  be  taken  into  account. 

v  1 .  =  -  ?qj  E} 

m  .  co , 

i  1 

... (3  .9) 

—2  i  = 

m  .  to  „ 

i  2 


vie 


1  9e  E  1  , 


mew  1 


V  o 
—2  e 


'9e  E2>; 


m  to0 
e  2 
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3 • 4  Derivation  of  Second  Order  Currents 


The  equations  for  the  second  order  currents  can 

be  derived  in  a  manner  similar  to  the  procedure  used  to 

obtain  equation  (2.12).  The  magnetic  field  Bq  has  no  effect 

on  either  P  o r  Q  (equations  (2.14)  and  (2.15)),  however, 

since  k  =  0,  §.  =  0  .  The  effect  of  B  is  an  additional 

o 

second  order  term  v „  xu  B  in  the  linearized  form  of  equa¬ 
te  -z  o  H 

tion  (3.6)  and  v.3  j  x  u_2  6  0  in  equation  (3.7).  Hence,  the 
second  order  equation  of  motion  for  electrons  is 


i  (  W "  i  ;  )  V ,  “  V  .  V 

e  i  '  — 3  e  e  i  — 3 


He  E  3  +  u  ,  v  0 

- —  <j>  3  e  r  e 

me 


u  v  -  , 

—  r  3  e  <j)  e 


+  /He  Pu 


,m, 


..(3.10) 


and  for  ions 


(to-iv  •  )  v  ,  •  -  v  •  v  0  =  H  ;  e  3  +  u  ,  v  0  .  fi  .  -  u  v  ,  .  . 

v  e  i  '  — 3  i  e  i  — 3  e  —  —  —  <J>  3  i  r  i  — r  3  i  <j>  1 

m  . 
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+  I  H  i  \  ^  P  u 

— r 


0  0  0 


(3.M) 


where 


qa^3axizBo) 


mac 


(u_r  v 


3  a  <J>  -<pv  3  a 


r  )  H  g  0  i 

V" 


iL(f)V3ar^a  ~  ELrv3a(f)^a 


P  is  defined  by  equation  (2.14). 
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Equations  ( 3  »  10)  and  ( 3  »  1  1  )  can  be  reduced  in  terms 
of  their  respective  components  to  yield  two  matrix  equations. 


( CO  -  i  V  e  j  ) 
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(D 

_ 1 

v  3  e  z 

q  e  /  me 

"  ^  e  i 
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ii 
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...(3.12) 

and 

[M] [V]  =  [R]  ...(3.13) 
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(3,16) 
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From  equation  (3.12),  and  assumi 


ng  a!  > 


v 


e  i 


v, 


noqev3ez 


+  n0q5v 


u2 

4tt 


E  3 

z _ 

(  (jj  -  I  V  g  j  ) 


...(3.17) 


Using  Mason's  Rule,  the  solution  to  equation  (3.13)  can  be 
written  as 


V 


< 


(-1)K+,R.C 


K  1 


...(3.18) 


a.  i_  _  .  _  . 

where  \l  K  is  the  k  1  element  of  the  matrix  [  V  ]  ,  |  M  j  ,  is  the 

determinant  of  the  matrix  [  M  ]  ,  and  C  K  .  is  the  cofactor  of 

R  in  the  matrix  formed  by  replacing  the  column  of  [  M  ] 

by  [R],  Hence,  from  equation  (3.18),  and  since 

j  ~  =  n  q  _  v  0  +  n  q  .  v  0  .  , 

*-3  oNe-3e  oM  i  — 3  i 

n  o  |  e  |  ^ 

j  3  r  =  |m|  I  ( -  t  )lR ,.  ( C  ,  ,  +  C2l)  ...(3.19) 

j3*  =  — T^T”i  =  i  (  "  1  )R  1  (  C3  1  +  Cl)l)  ...(3.20) 

The  cofactors  in  equations  (3.19)  and  (3.20)  are 
now  summarized  in  their  exact  form. 


c  1  2 
C1  3 

C  i  4 


2  2 
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2  1 
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23 

24 


v  .[fi.fi  +  (  OJ  -  i  V  .  )  (  CO  -  i  V  .  )  1  +  V  .  V  2  • 
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Also 
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(3.23) 

C  42  = 

C  2  4 

C  3  3  = 

C  1  1 
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C34  ' 

C  1  2 

J 
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M  |  =  (oo2-fig)(co2-fi|) 
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..(3.25) 


where 


v  =  vei  +  vei 


fi 


H 


fi  ?  fi. 


(3.26) 

(3.27) 


If  values  for  R%  are  substituted  into  equations 
(3*19)  and  (3.20),  the  result  is 
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3.28) 
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where 
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Ml  H 


.  .  .  (3.30) 


M 


[wpe(C13+C23)  +  wpi(Cl4+C24)] 


^  oo2^  1  +  1" 
I  M  | 


.  .  (3.31  ) 


M 


(  C  i  l  +  C  2  i  )  -  ^(C12  +  C22) 

me  11  ^  1  m  |  1 ^  L L 


i  ^  O)2  (  CU-  I  V  )  I  +  1  ‘ 

me  |  M  | 


..  (3.32) 
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(3.34) 
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^e?  1 


(3.35) 


If  the  real  parameters,  £]>  c2,  3),  32?  ?]»  ?2 

n  ]  and  n  2  »  are  defined  by 


a  |  -  C  |  +  i  3  |  > 


b  |  =  £ |  +  in 
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Then 


£l 

CO 


.2  ( n  2 

^  p  e  l  n  H 


.2) 


(co2-ft2)  (oo2-ft?) 


i2 

00 


O)2  cofi 

p  e  e 


(  CO  2  "  ft  p  )  ( oo 2  -  ft  ?  ) 


00 


CO 


np  [  (  002  -fi2  )  (co2  "ftf)  -  2  (  ^  §  -  00  2  )  2  ] 


pe 


cov 


( oo 2  -0,2  )  2  (  w2  -fj?  )  2 


-2 

00 

£i 

CO 


2w2e(fi^-oo2)cofiel  cov 

/  2  n2x  2  /  2  n2  \  2 

\  CO  "ftg  j  1  CO  -ft  |  ) 


]_  (X)peLQ2  1  +  1  ~ 

me  (co2  -ft2  )  ( oo2  -ft?  ) 


(3.38) 


£2 

oo 


1  ft  e  oo^2  e  [  co 2  (  1  +  m  3  )  -  ft  2  1  +  ] 


mew 


(oo2  -ft2  )  (co2  -ft?  ) 


00 


i  2  /  4  2  n  2  ,  -  2  n  ^  \  i  +  i  - 

1  OOpglcO  +C0  ft  g  1  “  ft  ^  J  I  I  cov 

rrig  ( co2 -ft2  )  2  ( co2 -ft?  )  2 


n2 

CO 


1  ftg^pe  [co^4  (  (2  +  mr)  1~  +  2m'j)-co2ft^(1+m2)  1++ft 


2  ^  ,  +  ^202-.+ 

i  JiH  1 


]  cov 


me  co 


(oo2-ft?)  2  ( oo2  -  ft?  )  2 


where 


CO 


p  a 


2  l  + 

00 

p  a 


(3.39) 


35 


3.5  Differential  Equations  For  Second  Order  Fields 

The  differential  equations  for  the  second  order 
fields  are  the  same  as  equations  (2.22)  with  k  =  0. 
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(3.42) 


From  equation  (3.17)  it  can  be  seen  that  j  ^  is 
proportional  to  E3  only.  Hence,  the  only  solution  to  equa 
tion  (3.42)  is  the  homogenous  solution,  E3  5  0.  By 
substituting  equations  (3.28)  into  equations  (3.40)  and 
(3.41),  the  set  of  equations  to  be  solved  becomes 
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It  may  be  noted,  E3  exists  as  a  result  of  coupling 
between  field  components  caused  by  the  magnetic  field  B0„ 
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3.6  Solution  of  Second  Order  Field  Equations 

The  following  analysis  is  similar  to  Art.  2.6.  A 
F ou r i e r - B e s s e 1  expansion  for  E|  consistent  with  the  boundary 
condition,  E|(r=rQ)  =  0,  is  used.  From  equation  (3  .  H)  ,  an 
expansion  of  first  order  Bessel  functions  appears  appropri¬ 
ate,  hence  it  is  assumed  that 


E|(r) 
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AKJ,  (rAj 


(3.45) 
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J  ,  (rQXK) 
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o  r 


1  K 
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t  h 

where  s,  is  the  k  zero  of  a  first  order  Bessel  function 

1  K 

Following  the  methods  of  Chapter  Two,  it  may 
readily  be  shown  that 
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The  solution  for  the  second  order  fields  is 
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Equations  (3-51  )  and  (3  -52)  are  now  written 
form  more  adaptable  to  the  following  sections. 
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f(r)  =  r  3  [  J i  (  k  ]  t  r ) J  i  (  k 1 2  r ) ] 
3  r 


...(3.58) 


3.7  Energy  Absorbed  by  the  Ions 


The  rate  of  energy  absorption  by  the 


on  s  ,  W  .  ,  is 


given  by 
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Reft*  E  3  +  T  *  .  .  E  2  ') 

^  J  3 i r  r  J  3 i b  T  ; 


(3.59) 


where  j  .  and  j  .  are  the  ion  components  of  the  total 
3  i  r  3  i  <t> 

current  with  source  terms  excluded.  Equations  analogous  to 
equations  (3.28)  can  be  written  for  these  components. 
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Combining  equations  (3.59)  and  (3.60)  gives 
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.  .  .  (3.63) 


Under  the  assumptions  of  equation  (3.37) 
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3.7.1  Optimization  of  Energy  Absorbed  by  the  Ions 

The  procedure  used  in  this  section  will  be  similar 
to  that  used  to  optimize  the  fields  in  Chapter  Two,  however, 
some  simplifying  assumptions  will  be  made. 

Since  both  E3  and  E3  are  proportional  to  1 /A  , 

<p  r  k 

2 

then  the  energy  dissipated  will  be  proportional  to  (  1  / A ^ )  . 

If  the  real  part  of  1/A^  can  be  set  to  zero  for  some  < ,  say 
k  =  1,  then  (  1 / A  ^ )  will  be  proportional  to  1 /v^ j  ,  hence  W. 
will  be  proportional  to  1/v  .  for  the  term  <  =  1.  Depend¬ 

ing  on  how  small  vgj  is,  a  considerable  peaking  of  the 
energy  spectrum  may  be  expected  as  Re(A^)  0.  In  contrast 
to  the  unmagnetized  plasma,  the  portion  of  A(  contributing 
this  effect,  referred  to  as  P  ,  is  a  function  of  1. 

"  (Hi.)2  1*2]  2  ...(3.66) 
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The  equivalent  term  for  the  unmagnetized  plasma,  since  it 
was  not  a  function  of  \,  could  be  factored  as  common  to  all 
terms  in  the  series  expansions  for  the  fields.  Thus,  under 
optimum  conditions,  all  terms  in  the  series  were  optimized. 
For  the  case  now  being  considered,  only  a  single  term  can  be 
optimized.  This  is  not  a  very  serious  characteristic  since 
the  convergence  of  the  series  is  rapid  and  higher  order 
terms  have  very  little  contribution. 

When  an  exact  expression  for  Wj  is  considered, 
the  result  does  not  readily  lend  itself  to  analysis,  hence 
it  will  be  assumed  that  v  .  is  sufficiently  small  that 
under  optimum  conditions,  all  terms  <  ^  \  ,  can  be  considered 

negligible  compared  to  the  term  <  =  \.  Equations  (3-5^)  and 

(3-55)  then  reduce  to 
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Consider  equation  (3*66).  Under  the  assumptions 


of  equations  (3-37): 
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. . . (3.70) 


If  the  optimizing  condition,  Refp^  =  0,  is  applied,  the 
following  dispersion  relation  results. 
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(  00  2  -n|  )  (  002  -  ft2  )  +  W2e(fi^-  00  2  ) 


. . .  (3.71) 


This  relation  establishes  an  upper  bound  on  oo .  As 
[  ( oo2 -ft2 )  ( oo2 -ft2  )  +  k'pe^H”^)]  +  0  ,  X  x  ->  »  .  But  Xx  must  also 
satisfy  equation  (3-47).  Assuming  oo  j  2  >>  w>  an<^  using 
equation  (3.3): 
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where 
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Hence  it  is  being  implied  that  rQ  ->■  0  and  00^  00 .  If  this 

limiting  condition  is  not  to  be  reached,  then 
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Substituting  equation  (3-71)  into  equation  (3 


gives 
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It  is  also  necessary  to  evaluate  the  term  ^  to 
since  A.  =  ( w  2  -  X2 j  P  x  . 
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T  ,  defined  by  equation  (3*53),  may  be  approx 
by  its  real  part. 
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This  result  is  consistent  with  equations  (3-37). 


Substituting  equations  (3.64)  and  (3.65)  into 
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equation  (3.69)  gives: 
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and,  using  equation  (3*78),  this  becomes 
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Consider  now,  the  expression  for  B,  (equation 

(3.56)).  B  is  proportional  to  bj.  [  l  +  aj  ]  +  £2^.1*  however, 
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this  term  can  be  approximated  by  its  real  part. 
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In  the  last  approximation  terms  'v  0(mr)  were  neglected, 
however,  at  this  stage  of  the  analysis,  the  approximation 
has  no  serious  effect  on  the  final  result. 

Hence,  using  equations  (3.56),  (3.75),  (3-76), 

(3.80)  and  (3.81),  the  rate  of  energy  absorption  by  the  ions 


becomes : 
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To  reduce  Wj  to  equation  (3.82)  from  equation  (3.80),  approx¬ 
imations  as  high  as  'v  0(4mr)  were  made  on  component  terms, 
however,  the  cummulative  effects  of  all  approximations 
tended  to  cancel.  In  addition,  none  of  the  approximations 
were  of  a  critical  nature. 


3.7*2  Maximization  of  the  Optimum  Solution 


Equation  (3.82)  may  be  considered  an  optimum 
solution  for  the  rate  of  energy  absorption  by  the  ions. 

The  optimizing  condition  prevailing  is,  Re(Pl)  =  0,  and 
this  specifies  relations  between  the  variables  r  ,  w,  w j  , 

B  and  nQ  (see  equations  (3.71)  and  (3.72)).  Therefore,  to 
maximize  W - ,  the  parameters  in  equation  (3*82)  must  be 
varied  in  accordance  with  the  dispersion  relation  and  within 
the  bounds  of  validity  of  the  analysis  -  i  .  e  .  :  co  ^  ^ 

sufficiently  large;  high  phase  velocities;  equations  (3*37) 
satisfied. 

Since  w  |  2  00  as  w  »  then  Wj  -*  0.  Now 

consider  gj  <  15ftj.  Then 
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Hence 


...(3.85) 
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W.  will  be  a  maximum  when  ft.  is  a  maximum  and 

i  i 

co  |  2  •  s  a  minimum.  If  shielding  of  the  incident  waves  by 

electrons  is  not  to  be  effective,  then  the  minimum  value  of 
co  |  2  '  s  roughly 
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The  dispersion  relation,  equation  (3-71),  is  satisfied 
when 


co 


^mn10 


...(3.88) 


Two  cases  will  now  be  considered.  In  the  first 
case  it  will  be  assumed  that  cope  >  fte  ,  and  in  the  second  it 
will  be  assumed  that  cope  <  fte. 
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If  this  is  to  be  a  valid  maximum,  then  there  is  an 
bound  on  o^mn  approximately  given  by 
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Otherwise,  equation  (3*86)  does  not  hold  since  go  > 
Assuming  equation  (3-90)  is  satisfied,  then 
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It  must  now  be  determined  which  term  in 

expression  for  the  energy  absorption  to  optimize  ( 

and  which  modes  of  excitation  should  be  used  (ie: 

It  is  desirable  to  make  F  as  large  as  is  practi 
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In  Appendix  2,  it  is  shown  that  for  a  given  m  and 

n,  Fimn  must  eventually  begin  to  decrease  and  approaches 

zero  as  i  approaches  infinity.  The  convergence  is  shown  to 

be  rapid  and  the  best  choice  for  i  is  i  =  1.  If  i  =  1, 

equation  (3-90)  is  satisfied  for  all  m  and  n,  m  ^  n. 

The  best  selection  for  m  and  n  is  governed  by 

practical  limitations.  The  tables  in  Appendix  2  indicate 

that  as  m  and  n  increase  for  a  fixed  value  of  i,  F 

*  \  mn 

increases,  however,  it  also  becomes  more  difficult  to  excite 
the  modes.  In  addition,  the  radius  of  the  waveguide  be¬ 
comes  large  as  m  and  n  increase,  since  when  W.  is  a  max¬ 
imum: 
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Figure  3*2  Illustrates  the  variation  of  F  for 

i  mn 

various  value  of  i,  m  and  n.  A  trend  is  evident  and  this 
is  indicated  by  continuous  lines  although  the  actual  points 
are  discrete.  Beyond  the  range  considered,  the  results  are 
practically  of  little  interest,  but  the  points  appear  to 
be  approaching  a  limiting  value  as  they  must  since  there  is 
only  a  finite  amount  of  energy  in  the  incident  waves. 
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Figure  3*2  Variation  of  F 
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magnetic  field  intensity  is  increased. 
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To  satisfy  the  inequality 
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If  inequality  (3*95)  holds,  then  equation  (3*86)  becomes 
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would  have  to  be  of  order  greater  than  25.  In  addition, 
the  radius  of  the  waveguide  would  be  excessive,  therefore 
this  case  is  not  a  practical  consideration. 

Hence,  case  [1]  is  preferable  to  case  [2]  and  the 
fact  that  a  lower  magnetic  field  is  used  is  an  added  practical 
advantage . 


3.7.3  Sensitivity  of  Wj  Near  the  Optimum  Point 


In  this  section  the  sensitivity  of  the  optimum 

solution  to  small  fluctuations  in  00,  n  and  B  will  be 

o  o 

investigated.  As  in  Art.  2.8,  this  will  be  done  by  consid¬ 
ering  the  term  from  which  the  optimizing  criteria  was 
derived;  in  this  case  P  ^ .  Defining  the  real  quantities 
and  Y  ^  by 


P 


x 


i  Y 
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then  an  optimum  condition  prevails  provided  Xx  does  not 
become  large  compared  to  Y;.  The  following  inequalities 
should  be  satisfied 
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To  simplify  the  analysis  it  will  be  assumed  that 
equation  (3.90)  holds  so  that  the  assumption,  oo  <  1  5  SI  j  ,  is 
valid.  Then,  from  equation  (3.75) 
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Defining  a  quality  factor,  Q ,  by 
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then 
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Therefore  variations  in  nQ  are  restricted  within 
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Therefore  variations  in  B0  are  restricted  within 
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3  .  7  .  Energy  Dissipation  in  a  Hot  Diffuse  Plasma 


In  order  to  give  a  clearer  indication  of  some  of 
the  actual  numbers  involved  in  ion  heating,  a  hot  diffuse 
plasma  will  now  be  considered.  It  will  be  assumed  that  the 
incident  waves  are  TEqj  and  TEQ2fnodes.  The  plasma  is 
characterized  by  the  following  parameters. 
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first  term  in  the  series  expansion  for  W.  wil 
hence  i  =  1.  It  will  also  be  assumed  that 
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Using  the  formulas  developed  in  Art.  3-7.2  for  the  case 
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W,  *  0  .  1  07F,  2  |  |  E  Je|  I  2  *  5.1xlO-'t|ElE||2 

rQ  s  7.24s,  ,  s  32.5  cm 

0  1  2  1 

B0  -  3.19  kilogauss 

ft  j  -  3 . 02  x 1 07  sec  ' 

a)  -  9.54x10^  ai21  ~  8.11x)07  sec"' 


From  Art.  3-7.3 
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An  instructive  way  to  consider  the  problem  is  to 
determine  the  approximate  time,  t,  the  sources  must  be 
turned  on  to  achieve  a  1%  increase  in  temperature.  Equating 
the  change  in  energy  of  the  ions  to  the  energy  absorbed  by 
them  in  time  t 


nomi ve i ( Ave i )  *  wiT 


...(3.112) 


The  power  output  of  each  source  is  given  by 
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The  integration  is  carried  out  over  the  cross-section  of  the 
waveguide  and  l  denotes  either  source  1  or  source  2.  From 
the  equations  of  Art.  1.3 


U 


s 

i 


...(3.114) 


Hence,  for  the  case  now  being  considered 

T  =  /  cro\2rtomi  ve  i  <Ayei>  jq(si  1 )  -J  0  ( s  |  2 ) 
\32  J  5 .  I  « 1  O'11  U^U| 


...(3.115) 


A  1%  increase  in  ion  temperature  requires  a  0.5%  increase 
in  the  thermal  velocity  of  the  ions,  therefore 

4  - 1 

A vQ  .  -  6.5*10  cm-sec 

U  I 

Then,  assuming  the  sources  have  equal  power  outputs 

t  *  3.25*1 026 

( Us  )  2 


To  achieve  a  1%  increase  in  temperature  in  1  second  under 
ideal  conditions,  sources  with  an  output  power  of  about 
1 8 0 0  kilowatts  would  be  required.  Since  the  output  of 
these  sources  would  have  to  be  pulsed  the  time  will  be 
correspondingly  increased.  It  will  be  further  increased  if 
the  required  value  of  Q  cannot  be  realized.  Hence,  the 
problem  as  treated  here,  under  the  assumptions  made,  yields 
a  low  rate  of  heating,  however  it  should  be  possible  to 
improve  these  results  by  several  orders  of  magnitude  and 
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suggestions  to  this  effect  follow  in  Art.  3.8. 

As  a  final  point  of  interest,  it  will  now  be 
checked  that  the  high  phase  velocity  approximation  is  valid. 
Since  k  =  0,  the  phase  velocity  in  the  z-direction  is 
infinity  hence  the  condition  (equation  (1.4))  is  certainly 
satisfied  in  this  direction.  Now  consider  the  direction 
perpendicular  to  z. 

X 

i 

Hence 

co  ~  S7  j  q  7 

- - - L  4 . 3  5  x  1  O'*  >>  Vqj  -  1  .  3  x  1  0  ' 


12 

e  1  a  5  .  2x 1 0  1  z 


X 


6  e 


5.5x10 


8 


Therefore,  the  high  phase  velocity  approximation  is  justi¬ 
fied  for  this  case. 


3.8  Concluding  Remarks 

Although  the  degree  of  heating  that  can  be  achieved 
is  small  for  the  problem  as  considered  here,  this  should  not 
be  accepted  as  indicative  of  the  overall  potential.  A 
principal  factor  contributing  to  the  low  value  of  energy 
absorption  is  the  large  radius  that  is  required  to  allow  a 
cut-off  signal  (k  =  0),  which  leads  to  the  consequent  high 
power  output  of  the  sources  that  is  necessary  to  achieve 
appreciable  heating.  If  this  condition  is  relaxed,  then 
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the  radius  can  be  set  arbitrarily.  Although  the  available 
energy  is  not  being  used  as  efficiently  as  possible,  the 
amount  of  heating  that  is  realizable  will  be  greater  since 
a  smaller  confining  structure  can  be  used.  Thus  for  sources 
of  a  given  power  output,  the  energy  density  in  the  structure 
will  be  greater.  A  possible  method  to  decrease  the  loss 
in  efficiency  for  cases  where  k  ^  0  would  be  to  employ  a 
toroidal  or  figure  eight  (  s t e  1  1  a r a t o r )  structure.  A  further 
advantage  with  smaller  radii  is,  the  area  across  which  the 
magnetic  field  must  be  applied  is  less. 

The  reason  for  selecting  cut-off  signals  was  based 
on  the  results  of  an  analysis  carried  out  (6)  for  an 
infinite  plasma  excited  by  plane  waves.  The  wavevector  of 
the  second  order  fields  was  defined  by  a  component  perpen¬ 
dicular  to  the  direction  of  the  static  magnetic  field,  k^  , 
and  a  parallel  component,  kM  .  It  was  found  that  setting 
k(l  to  zero  gave  the  better  results.  It  is  interesting  to 
compare  this  problem  to  the  problem  considered  here.  The 
ratio  of  the  rate  of  energy  absorption  in  the  infinite 
plasma,  W? ,  to  the  value  of  Wj  given  by  equation  (3*86)  is: 
W?/W.  =  A/F  .  In  general,  this  ratio  is  large,  as  can  be 

seen  from  the  values  of  F  given  in  Appendix  2.  This 

i  m  n  3  r  r 

difference  can  be  related  to  the  equations  for  the  second 


order  fields  and  arises  because 
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/  r  J  |  (  r  A  ) 
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3  r  r  3  r 


2p" 


d  r 


r°  2 

/  r  J  |  (  r  A  t  )  d  r 
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x  J  1  (  r  X 


<  < 


where  ( 6 ) 


=  -  i  I  e  I  kj_  (  E  1  •  E  2  * ) 


<*>  1  to  2 


In  an  Infinite  plasma  excited  by  plane  waves,  the  degree  of 
coupling  between  the  incident  waves  is  constant  throughout 
and  P*5  has  no  spatial  dependence.  However,  in  a  finite 
plasma  of  the  type  considered  here,  P  does  have  a  radial 
dependence  which  adversely  affects  the  net  amount  of  coup¬ 
ling.  A  possible  explanation  is  as  follows: 


Consider  each  incident  wave  as  imparting  a  component  of 
electromagnetic  momentum  to  the  second  order  wave.  For  the 
momentum  transfer  to  be  a  maximum,  these  components  must 
add  over  all  space;  the  resultant  second  order  field  being 
proportional  to  the  spatial  average  of  all  momentum  contri¬ 
butions.  In  the  infinite  plasma,  these  components  are  fixed 
and  give  the  same  contribution  over  all  space,  however,  in 
the  finite  plasma,  these  components  vary  in  magnitude  and 
direction  and  the  net  effect  upon  averaging  is  small. 
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The  underlying  reason  for  the  difference  in  the 
degree  of  coupling  between  the  finite  and  infinite  case  is 
in  the  nature  of  the  waves  which  excite  the  plasma.  In  the 
infinite  plasma  we  are  concerned  with  the  coupling  of 
travelling  waves  in  the  direction  perpendicular  to  the 
magnetic  field,  whereas  in  the  finite  plasma,  we  are  coup¬ 
ling  standing  waves.  The  existence  of  standing  waves 
specifies  a  discrete  (since  the  modes  of  excitation  can 
be  adjusted  discretely)  phase  relation  between  the  incident 
waves.  To  achieve  travelling  wave  coupling  in  a  finite 
structure,  the  magnetic  field  would  have  to  be  applied 
across,  rather  than  along,  the  structure.  Then  the  incident 
waves  across  the  magnetic  field  appear  as  travelling  waves 
as  in  the  infinite  plasma.  This  particular  point  is  one 
which  deserves  detailed  investigation.  A  further  area  open 
to  study  is  co n s i d e r a t  i  on  of  other  types  of  geometry  for 
the  confined  region. 

The  sensitivity  of  Wj  to  small  fluctuations  in 
any  of  the  parameters  is,  at  present,  somewhat  restrictive. 
To  control  the  density  of  the  plasma  to  one  part  in  10^  is 
very  difficult  and  is  a  further  justification  for  the  use 
of  smaller  structures.  The  quality  factor,  Q,  is  also  high. 
The  hot  diffuse  plasma  considered  in  Art.  3-7-^  required  a 
Q  of  *4 . 6 x  1  0  7  t  however,  a  more  desirable  figure  would  be 
about  105  since  this  is  more  within  the  realm  of  present 


technology. 
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It  was  mentioned  in  Art,  3 .7 » 4  that,  since  the 
power  levels  involved  to  achieve  significant  heating  were 
high,  the  sources  would  have  to  be  pulsed  rather  than 
having  a  continuous  output.  This  disadvantage  could  be 
slightly  lessened  by  using  several,  rather  than  just  two, 
sources. 

Assuming  a  method  similar  to  the  type  being 
considered  was  to  be  used  for  ion  heating,  it  should  be 
noted  that  the  rate  of  heating  increases  as  temperature 
increases  since  the  collision  frequency,  v  •  ,  decreases 
(or  at  worst  the  rate  remains  constant,  since  sensitivity 
increases  in  direct  proportion  to  1  / v  e  j )  ,  This  is  in 
contrast  to  the  method  of  i o n - cy c 1 o t r o n  heating  discussed 
by  Berger  et  al  (9),  where  at  high  temperatures,  the  rate 
of  heating  decreased  as  temperature  increased. 


. 


CHAPTER  A 


DISCUSSION 


The  analysis  presented  in  this  thesis  indicates 
a  suitable  approach  to  some  aspects  of  the  problem  of 
nonlinear  interaction  of  high  frequency  waves  in  bounded 
plasmas,  once  basic  assumptions  as  to  the  nature  of  the 
plasma  and  what  is  to  be  done  with  it  are  made.  In  theory, 
the  problem  could  be  generalized  to  a  greater  extent  by 
assuming  an  arbitrarily  shaped  boundary  and  expanding  the 
fields  in  some  set  of  orthoganal  functions,  $  ( ip  K  ,  r_)  . 

(In  this  thesis  $(i|^K,r_)  ->  Jn(XKr)).  The  main  problem  in 
practice,  of  course,  would  be  to  find  a  tabulated  set  of 
orthoganal  functions. 

One  aspect  of  the  problem  which  requires  research, 
and  to  the  author's  knowledge  has  not  yet  been  investigated 
is  wave  coupling  in  a  bounded,  non-uniform  plasma.  This  is 
a  more  realistic  model  and  it  is  well  known  that  such  a 
plasma  will  in  general  support  more  than  one  normal  mode  of 
oscillation  (10),  (11),  whereas  the  uniform  plasma  has  only 

one  normal  mode  at  the  plasma  frequency.  Another  factor 
which  contributes  to  additional  resonances  is  the  thermal 
velocities  of  the  electrons  and  ions  which  must  be  taken 
into  account  if  phase  velocities  are  small.  A  further 
complicating  factor  which  now  arises  due  to  these  additional 
considerations  is  that  of  correlating  the  field  spectrum  to 
the  density  distribution  if  wave  mixing  is  to  be  used  as  a 
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diagnostic.  This  is  a  point  where  the  mixing  of  laser  beams 
enjoys  an  advantage  over  the  general  mixing  of  waves  through¬ 
out  the  plasma.  With  lasers,  the  interaction  can  be  inves¬ 
tigated  in  very  local  regions  across  which  density  distri¬ 
bution  can  be  ignored. 

With  respect  to  the  section  on  ion  heating, 
although  the  rate  of  energy  absorption  by  the  ions  was  small 
for  the  case  considered,  the  investigation  has  indicated 
the  areas  to  which  future  studies  should  be  directed.  A 
major  limiting  factor  for  the  amount  of  heating  that  can 
be  achieved  is  the  necessarily  high  frequency  of  the 
incident  waves  to  permit  them  to  thoroughly  penetrate  the 
plasma.  The  amount  of  heating  varies  inversely  as 
Not  only  would  the  heating  increase  substantially  if  lower 
frequencies  could  be  used,  but  a  considerable  reduction  in 
Q  may  be  expected.  In  addition,  the  method  can  then  be 
applied  to  higher  density  plasmas  with  direct  application 
to  thermonuclear  heating.  The  problem  of  using  lower 
frequency  waves  is  currently  being  investigated  by  C.  E. 
Capjack  (12). 

In  conclusion,  there  are  yet  many  aspects  of  the 
problem  of  nonlinear  coupling  of  waves  in  plasmas  that 
require  detailed  study.  Two  major  applications  have  been 
mentioned  and  there  are  perhaps  many  more.  In  particular, 
the  achieving  of  significant  amounts  of  ion  heating  has 


great  potential. 
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APPENDIX  1 


DERIVATION  OF  DIFFERENTIAL  EQUATIONS  FOR 


SECOND  ORDER  FIELDS 


In  this  Appendix,  differential  equations  for  the 
the  second  order  fields  of  the  unmagnetized  plasma  (Chapter 
2)  are  derived.  The  results  are  directly  applicable  to  the 
magnetized  plasma  (Chapter  3). 


From  equations  (1.5)  and  (1.8) 


v*i_3  =  “  iL°p3 

and,  from  equations  (1.7)  and  (1.8) 
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Substituting  equation  ( A 1  .  1 )  into  equation  (1.5)  gives 


Hence 
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The  components  of  V  E J  are  given  by 
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Substituting  equations  (A1.4)  and  (A1.5)  into  equations 
(A1.2)  results  in  equations  (2.2^).  If  k  is  set  to  zero, 
then  equations  (3.^0)  -  (3.^2)  result. 
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THE  FUNCTION  F; 


mn 


i  mn 


1  T 


£- 


i  1 


+  m  -  S  1  ;  I 


r5 1  i  1  i  2 


2  J ) ( rXj ) 
Jo ( ro A  i ) 


. . . (A2 . 1 ) 


It  will  now  be  shown  that,  1 i m  F j mn  =  0,  for  any 
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Since  G(r)  and  H(r)  both  satisfy  Dirichlet's 
con d  i  t  i  on s ^ i n  the  closed  interval  [0,rQ],  for  any  finite 
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value  of  m  and  n  (recall  that  k  ^  j  =  s]m/r0>  and  ^ )  2  =  sln^ro^’ 
then 

I  °°  J  |  (  r  X  .  )  ro 

7 2  l  V  /rJ|  (rA;)G(r)dr,  and 

roi=l  ^oro^i  0 
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— j  I  ~2~, - —  /  r  J  |  (  r  X  j  )  H  (  r  )  d  r  ,  are  convergent  ( 1  3 )  . 

roi=l  ^  0  (  roM  '  0 

Hence,  by  comparison,  the  series  in  equations  (A2.3)  and 
(A2.5)  are  convergent.  Therefore 

00 

l  [  F  i  mn  -I  is  bounded,  =>  lim  F  ;  mn  =  0 
i  =  1  i 


A  check  of  Tables  A2 . 5  ■  A2.8  will  verify  that  the  conver¬ 
gence  is  rapid  and  that  the  best  choice  for  i  such  that 

F.  is  maximized  is  i  =  1. 

1  mn 

2 

Values  of  I  .  ,  and  m  s , . I . «  are  given  in  tables 
1  I  r  I  1  1  2  3 

A2.1  -  A2.4.  The  integrals  were  solved  numerically  on  an 

IBM  70A0  computer  by  dividing  the  integrands  into  fifty 
segments  and  applying  Simpson's  rule.  Those  numbers  which 
are  of  questionable  accuracy  due  to  computer  roundoff  error 
and  error  inherent  in  Simpson's  rule  are  separated  from  the 
other  numbers.  The  integrals  are  given  by 
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To  facilitate  a  computer  sol ut ion  f  the  derivative  terms  were 
expanded . 
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Note  that  equations  (A2.9)  and  (A2.10)  verify  that  G  (  r )  and 
H ( r )  satisfy  Dirichlet's  conditions. 


A  function  f(x)  satisfies  Dirichlet's  conditions  in  the 
interval  [a,b]  if  (13): 

1.  f(x)  has  only  a  finite  number  of  maxima  and  minima  in 
[  a  ,  b  ]  . 

2.  f(x)  has  only  a  finite  number  of  finite  discontinuities 

in  [a,b]  -  and  no  infinite  discontinuities. 
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